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We report on a search for the super symmetric partner of the top quark (stop) produced in tt events 
using 110 pb -1 of pp collisions at yfs = 1.8 TeV recorded with the Collider Detector at Fermilab. In 
the case of a light stop squark, the decay of the top quark into stop plus the lightest supersymmetric 
particle (LSP) could have a significant branching ratio. The observed events are consistent with 
Standard Model tt production and decay. Hence, we set limits on the branching ratio of the top 
quark decaying into stop plus LSP, excluding branching ratios above 45% for a LSP mass up to 40 
GeV/c 2 . 
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With the observation in 1995 of a heavy top quark [J, , an important prerequisite was met 
for low energy supersymmetry (SUSY) || to explain electroweak symmetry breaking. In the 
Minimal Supersymmetric extension of the Standard Model (MSSM), all known particles of the 
Standard Model (SM) acquire supersymmetric partners, or superpartners. For fermions the 
superpartners are bosons, while for bosons the superpartners are fermions. We assume con- 
servation of a multiplicative quantum number 7£-parity, which requires these new particles to 
be produced in pairs and prevents decays of the lightest supersymmetric particle (LSP). From 
cosmological considerations 0, the LSP is normally assumed to be the lightest neutralino. 

The large Yukawa coupling of the top quark opens up the possibility of a large mass 
splitting in the third generation of fermionic superpartners (the squarks and sleptons). The 
superpartners of the right-handed and left-handed top quark (the stop squarks) combine to 
form the mass eigenstates. The lightest stop squark (ti) could then be lighter than the super- 
partners of all other squarks. Most limits on squark masses || do not apply to the stop squark 
because they are usually based on a model of five degenerate squarks. Current stop squark 
mass limits are significantly lower |J than these limits or based on the assumption of a very 
heavy chargino (xf) @. The latter searches are complementary to the analysis presented here 
since the stop decay mode c + LSP does not coexist with the decay mode b + xf ■ If the stop 
squark is light, decays of the top quark into stop plus the lightest neutralino could be kinemat- 
ically allowed. If this neutralino is the LSP it will be stable and only weakly interacting. Such 
a particle would pass through the detector without interaction, causing a considerable energy 
imbalance. For the stop squark we assume decays analogous to the Standard Model top quark 
decay, i.e. into chargino and 6-quark. The chargino could then decay into a LSP plus either a 
quark-antiquark pair or a lepton and neutrino. 

Branching ratios as large as 40 to 50% for the top decay into stop have been suggested M. 
In such scenarios about one half of it events would have one SM and one supersymmetric top 
decay. If the SM top decay to a W and a 6-quark is followed by the leptonic decay of the W, 
then the selection criteria of the online leptonic trigger, the offline dataset selection and the 
it event identification |]J will all be satisfied. The decay of the second top quark can then be 
used to search for decays into stop. 

The CDF detector |IJ is well suited to search for supersymmetric top quark decays. The 
following components are relevant to this analysis: the central tracking chamber, which is inside 
a 1.4 T superconducting solenoidal magnet, measures the momentum of charged particles with 
a resolution of 5]?t/pt = 0.001 * p T (p T in GeV/c) [|TTJ. The silicon vertex detector, with 
an inner radius of 3 cm and an outer of 8 cm, identifies secondary vertices with a resolution 
of 130 jum in the transverse plane. The electromagnetic and hadronic calorimeters cover the 
pseudorapidity region \r)\ < 4.2 and are used to identify jets and electrons, and to measure 



the missing transverse energy ID- 11 - An outer layer of drift chambers provides muon 

identification in the region \rj\ < 1.0. 

The search reported here is based on llOpb" 1 of pp collisions at y/s = 1.8 TeV recorded 
with the Collider Detector at Fermilab during the 1992-93 and 1994-95 collider periods. The 
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analysis is a combination of a CDF single lepton plus jet top analysis and a kinematic 
analysis ||12|| . The analysis cuts are slightly revised to improve sensitivity 

The leptonic W decay from the SM top decay yields an energetic lepton. Events are 
selected as in the single lepton plus 6-jet top analysis by requiring a central electron (|?y| < 1.1) 
or central muon (\rj\ < 1.0) with transverse momentum > 20GeV/c. 

The neutrino from the W decay, as well as any LSP's, will escape the apparatus without 
detection, resulting in an energy imbalance. The i^x measured by the calorimeter is corrected 
if the lepton is a muon. We have increased the E^ requirement, from Et > 25 GeV in the 
single lepton plus 6-jet top analysis to E^ > 45 GeV, as our signal is expected to have a 
harder E^ spectrum. To reject non-W background, we require the transverse mass Mp of 
the lepton and E^ system to be larger than 40 GeV. While the top analyses are concerned 
with separating it from W plus jet production, in this analysis we are interested in minimizing 
W plus multijet background and then focusing on separating SUSY top decays from SM top 
decays. In both SUSY and SM top decays the W has a substantial transverse momentum. An 
additional requirement that the lepton— En? system have > 50GeV/c has therefore been 
made. 

In addition to the 6-jet from the SM top decay, we have three additional jets from the 



other top decay when the W or the chargino decays hadronically. We require two jets [13| with 
transverse energy > 20 GeV and a third jet with > 15 GeV all within \rj\ < 2.0. Jet 
energy is corrected according to an average response function prior to event selection [ 12 j . We 
require one of the jets to be identified as a 6-jet candidate. We use a secondary vertex tagging 
method [1], based on SVX information, that reconstructs secondary vertices from B hadron 
decays. As in the kinematic top analysis [pL2[| , we require the three jets to have large polar 
angles in the rest frame of the lepton, E T , and jets: | cos(^*(jet i ))| < 0.9, | cos(0*(jet,-))| < 
0.8, and | cos(0*(jet fe ))| < 0.7, where the jets i,j,k are ordered according to |cos(6 l *)|. These 
requirements were relaxed from the requirement of | cos(6*(iet i j k ))\ < 0.7 in [l2j in order to 
obtain a good acceptance for the SM-SUSY top decays. For simplicity, these polar angles 
are calculated assuming a null longitudinal component for the neutrino. In order to insure 
that the three jets are well separated in rj-cj) space, we require Ai?(jet, jet) > 0.9, where 
Ai?(jet, jet) = a/At/ 2 + A0 2 is the minimum distance between the jets. 

With these jet requirements our sample is defined. In the data 9 events pass the above 
cuts. For the theoretical it production cross section of 5pb JE4J we expect 9.5 events. Non-top 



Standard Model background contributes less than half an event. 

Jets from a SUSY top decay are less energetic than jets from a SM top decay due to the 
presence of LSP's and possible small chargino mass. In order to best distinguish between SM 
and SUSY top decays, we combine the transverse energy information of the second and third 

most energetic jet in a likelihood variable defined as: R L = vSM _ SUSYiE jet2 )xpSM _ SUSY J Ej et3y 

where V are the expected differential transverse energy distributions ^j§^ evaluated from the 
Monte Carlo. 
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Figure 1: Transverse Energy distribution of 
a) the second E^ 1 (2) and b) third jet E^ (3) 
for events with two SM top decays and for 
events with one SM and one SUSY top de- 
cay. 



Figure 2: a) Comparison o/ln(i?^) for SM- 
SM top decays and SM-SUSY top decays af- 
ter all the cuts have been applied, b) The 
ln(i?i) distribution for CDF Run 1 data. 



Figure 1 shows those Et distributions while Figure 2 shows the distribution of the likelihood 
variable for SM and SUSY top decays. Events with one top quark decaying into stop plus LSP 
are clustered in the region of negative values of ln(i?^), whereas events with two SM top 



decays are at positive values The region ln(i?x) < — 1 defines our SUSY search region. 
The region ln(i?^) > — 1 is dominated by double SM top decays and will be used to normalize 
the expected number of these decays. Our search is then independent of the it production 
cross section |L6 . 



To study the distributions of the kinematic variables in the super symmetric top decays, 
the signals are produced with the ISAJET Monte Carlo generator |L7| and passed through 



detector simulation programs [|18|]. We have fixed the top quark mass to 175 GeV/c 2 and varied 
stop, chargino, and LSP masses. Different branching ratios for top decaying into stop plus LSP 
are obtained by analyzing the three possible combinations of top decays (SM-SM, SM-SUSY, 
and SUSY-SUSY) individually and recombining them with appropriate weights. 

Systematic uncertainties in the simulation of it events can impact both the expected number 
of events and the shape of the two jet Et distributions and thus the likelihood distribution. 
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Figure 3: Top into stop plus LSP branching ratio limits as a function of stop and chargino 
masses for a LSP mass of 20 GeV/c 2 (left) and 40 GeV/c 2 (right). 



We have evaluated the systematic uncertainties for multiple points in the parameter space of 
the stop and chargino masses. The systematic uncertainties are expected to become significant 
close to the kinematic bounds. In the region of small chargino mass the uncertainties due to 
gluon radiation and the calorimeter energy scale are important. The procedure used to evaluate 
these systematic uncertainties is the same as that of Ref. [0. The uncertainty due to the it 
production cross section and to the integrated luminosity is negligible since we normalize the 
Monte Carlo predictions to the data in the SM dominated region of large likelihood. We note 
that this normalization has a large statistical uncertainty due to the small number of events in 
this region. The uncertainty in the mass of the top quark becomes important in the region of 
large stop mass, close to the kinematic limit. The effect of using parton distribution functions 
other than CTEQ-3 (LO) || is within the statistical uncertainty of the Monte Carlo samples. 
The total systematic uncertainty for both expected number of events and fraction of events 
with hi(R L ) < —1 from SM-SUSY top quark decays is typically around 40%. The major 
systematic uncertainty comes from the calorimeter energy scale and it varies between ±15% 
and ±25% in the region of the analyzed parameter space. The uncertainty on the top quark 
mass contributes about ±15% for most of the parameter space. The uncertainty due to the 
gluon radiation, on both number of events and shape, is always less than ±10%. 

All nine events observed in the data cluster in the SM-like region of large ln(i?^). To set 



8 



a limit on the branching ratio of top decaying into stop plus LSP, we calculate the branching 
ratios that would yield at least one event in the SUSY-like region 95% of the time as a function 
of stop, chargino and LSP mass. The method used is essentially a Bayesian-style integration 
over the systematic and statistical uncertainties in the SM-SM, SM-SUSY, and SUSY-SUSY 
contributions, where the uncertainties are assumed to be Gaussian distributed. Figure |3| shows 
the 95% confidence level top decaying into stop plus LSP branching ratio limit as a function 
of stop and chargino mass, for a LSP mass of 20 and 40GeV/c 2 . For larger LSP masses the 
kinematically allowed region shrinks. The sensitivity of this analysis, however, stays rather 
constant. 

In conclusion, we have looked for the top decay into stop plus the LSP. In the case of 
a light stop squark, the top is allowed to decay into stop plus a LSP. The number of events 
observed in the data is consistent with the Standard Model top decay expectation. We exclude 
branching ratios for top decaying into stop above 45% for an LSP mass up to 40 GeV/c . The 
upcoming run at the Tevatron, which will feature an approximately 20-fold increase in the 
total integrated luminosity as well as a significantly improved CDF detector, should allow 
these results to be greatly extended. 
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